Particle partitioning entanglement in itinerant many-particle systems 
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For itinerant fermionic and bosonic systems, we study 'particle entanglement', denned as the 
entanglement between two subsets of particles making up the system. We formulate the general 
structure of particle entanglement in many-fermion ground states, analogous to the 'area law' for 
the more usually studied entanglement between spatial regions. Basic properties of particle entan- 
glement are uncovered by considering relatively simple itinerant models. 



Introduction — In recent years, concepts from quan- 
tum information have proved useful for condensed mat- 
ter systems. One prominent example is the study of 
the entanglement between a part (A) and the rest (B) 
of the many-particle system, measured by the entan- 
glement entropy Sa- The entanglement entropy, Sa — 
— tr [pa In pa] , is defined in terms of the reduced density 
matrix pa = trg p obtained by tracing out B degrees of 
freedom. 

To define a bipartite entanglement, one has to first 
specify the partitioning of the system into A and B. 
The most commonly used scheme is to partition space, 
e.g., partition the lattice sites into A sites and B sites. 
However, for itinerant particles, with the wavefunction 
expressed in first-quantized form, one can meaningfully 
partition particles rather than space, and calculate en- 
tanglements between subsets of particles. Since each par- 
ticle has a label in first-quantized wavefunctions, indis- 
tinguishability does not preclude well-defined subsets of 
particles. Note that, with such partitioning, A or B do 
not correspond to connected regions of space. 

The distinction between particle entanglement and 
spatial entanglement was made relatively recently 0, 0, 
6j. In work reported since then, particle entanglement 
has been shown to be a promising novel measure of cor- 
relations 0, 0, 0, 0, 0] • In fractional quantum Hall states 
this type of entanglement reveals the exclusion statistics 
inherent in excitations over such states 0, 0|- Similar 
insight arises from particle entanglement calculations in 
the Calogero-Sutherland model [4j . For one-dimensional 
anyon states, particle entanglement is found to be sensi- 
tive to the anyon statistics parameter 0, 0] • 

Clearly, entanglement between particles in itinerant 
systems is a promising new concept, potentially useful 
for describing subtle correlations and the interplay be- 
tween statistics and interaction effects. A broad study of 
the concept and its utility is obviously necessary. Unfor- 
tunately, particle entanglement has till now been studied 
mostly in relatively exotic models, so that the literature 
lacks simple intuition about these quantities. This Letter 
fills this gap. We provide results for the simplest nontriv- 
ial itinerant fermionic and bosonic models, and present 
generic behaviors by generalizing available results. 



We first present upper and lower bounds for the en- 
tropy of entanglement S n between a subset of n parti- 
cles and the remaining N — n particles. We formulate 
a 'canonical' asymptotic form for fermionic and some 
bosonic systems. We next present results for a two-site 
Bose-Hubbard model. Through this toy model, we iden- 
tify two general mechanisms of obtaining nonzero particle 
entanglement in many-particle models. One mechanism 
is simply that of (anti-)symmetrization of wavefunctions, 
while the other is due to the formation of 'Schrodinger 
cat'-like states. The second mechanism is shown to be 
fragile, in the same sense that cat states are fragile in 
macroscopic settings. We next switch to true lattice mod- 
els, focusing on spinless fermions on a one-dimensional 
(ID) lattice with nearest-neighbor repulsion, sometimes 
known as the t-V model. We find similar mechanisms 
at work, in a nontrivial setting. In addition, our study 
of the t-V model enables us to present generic intuition 
about particle entanglement in many-particle systems, 
expressed in our canonical asymptotic language. 

Bounds - - A generic itinerant lattice system has 
N particles in L sites; we consider bosons or spinless 
fermions so that N < L. In every case, a natural upper 
bound for S n is provided by the (logarithm of the) size of 
the reduced density matrix pa = Pn, i.e., the dimensions 
of the reduced Hilbert space of the A partition. This size 

= C(L,n) for fermions and C{L — 1 + n, n) for 

bosons. The actual rank of p n can be much smaller due 
to physical reasons, so that the entanglement entropies 
are usually significantly smaller than the upper bounds, 
as we shall see in the examples we treat. 

In a bosonic system, S n can vanish, since a Bose con- 
densate wavefunction is simply a product state of indi- 
vidual boson wavefunctions, each identical. For fermions, 
however, anti-symmetrization requires the superposition 
of product states; for free fermions this causes p n to have 
C(N, n) equal eigenvalues. This provides a nonzero lower 
bound for S n in a fermionic system. 

Bosons : < S n < In C(L - 1 + n, n) , (1) 
Fermions : In C(N, n) < S n < In C(L, n) . (2) 

Canonical form — For large fermion number, N^>1, 
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we propose the following canonical form for the entan- 
glement of n<cTV fermions with the rest: 

S n {N) = lnC(N,n) + a n + 0(1/N) 
= n\nN + a' n + 0(1/N) . 

This form is suggested by results reported in Refs. [2, 0, 
, and in this Letter. For example, a n — nlnm for the 
Laughlin state at filling v = 1/m 2]. The same canon- 
ical behavior seems to hold for bosonic systems which 
lack macroscopic condensation into a single mode, e.g., 
bosonic Laughlin states Q , or hard-core repulsive bosons 
in one dimension Q. 

Subtle correlation and statistics effects can be con- 
tained in the behavior of the 0(1) term a n , and some- 
times also the 0(1/ N) term. Our calculations provide 
important intuition about how such effects show up in 
a n , as we summarize at the end of this Letter. 

Note that, for lattice sizes larger than TV, the generic 
behavior (|3|) indicates that the entanglement entropy 
does not saturate the upper bound (jT]) or (|2j) obtained 
from the size of the reduced Hilbert space. 

Two-site Bose-Hubbard model — We start with a toy 
lattice model, with only two sites. We will consider TV 
bosons on this 'lattice', subject to a Bose-Hubbard model 
Hamiltonian, to elucidate the basic mechanisms by which 
an itinerant quantum system can possess particle entan- 
glement. The Hamiltonian is 

H = -(blk+bik) + \U (blblhh + Sp|S 2 S 2 ) . 

(4) 

For U = 0, the system is a non-interacting Bose 
condensate, with each boson packed into the state 
(|1) + |2)). In the U — > +oo case, the system is a 
Mott insulator, with half the particles in site 1 and the 
other half in site 2. Such a state is simple in the "site" 
basis (second-quantized wavefunction) , but involves sym- 
metrization in the "particle" basis (first-quantized wave- 
function), leading to nonzero particle entanglement en- 
tropy. Finally, the U — > — oo limit involves all particles 
in either site 1 or site 2. The ground state is a linear 
combination of these two possibilities, which for large TV 
is a macroscopic 'Schrodinger cat' state. Such a state is 
somewhat artificial, because an infinitesimal energy im- 
balance between the two states will 'collapse' this state. 
For example, a 'symmetry-breaking' term of the form 
e b\b\, added to the Hamiltonian would favor site 2 
and destroy the cat state. The resulting state is a prod- 
uct state with zero particle entanglement. 

Incidentally, a 2-site model with off-site interaction V 
(instead of on-site U) has similar physics, with negative 
(positive) V playing the role of positive (negative) U. 

Two bosons in two sites — There is only one way of 
partitioning two particles (n = 1), so the only S n is Si. 
We expect Si = at U = 0, and maximal entanglement 
Si = In 2 for both 'Mott' state at U = +oo and the 
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FIG. 1: (Color online.) Entanglement between two bosons, 
in the ground state of a two-site lattice model with on-site 
repulsion. Solid curve is for the basic Bose-Hubbard model. 
Dashed (e = 0.1) and dash-dotted (e = 0.1U) curves illustrate 
fragility of 'cat' state via a eb[bi term. 

'Schrodinger cat' state at U = — oo. The Hilbert space is 
small; one can diagonalize the problem and calculate Si 
analytically as a function of U. We find Si(U) = Si(— U) 
interpolating smoothly between zero and In 2 ~ 0.6931 in 
both positive and negative directions (Fig.[TJ) 

We also demonstrate the fragility of the cat state by 
showing the effect of an eb\bi term. There is no appre- 
ciable effect for U > 0, but for U < the cat state is 
destroyed and we get Si — > oo for U — > — oo. 

Many bosons in two sites — For TV bosons, it is mean- 
ingful to study S n with n > 1. Labeling the basis 
states by site occupancies, i.e., as | iVi , A^a ) » the Mott 
and cat ground states are respectively \N/2, N/2) and 
{\0,N) + \N,0))/V2. The n-particle reduced Hilbert 
space has dimension n + 1; the reduced-space basis states 
can be labeled by the number of A bosons in site 1. 
In the Mott state \N/2,N/2), only the diagonal ele- 
ments of p n are nonzero and they are all equal; hence 
S n (U — * oo) = ln(n + 1). In the cat state, only 
two elements are nonzero, both on the diagonal; hence 
S n (U — * — oo) = In 2, independent of n. Fig. [2] demon- 
strates, via calculation from wavefunctions obtained by 
numerical diagonalization, that S n increases to ln(n + 1) 
and In 2 in the U — > ±oo limits. 

Both p n (U) and S n (U) can be understood in greater 
detail using approximations available in the literature 
9]. For U > 0, the coefficients of the ground 

state \GS) = ^N t \Ni, N — Ni ) can be approxi- 

mated by a gaussian ty Nl oc exp [(Ni - ±/V) 2 /ct 2 ], with 
No-~ 2 = (1 + UN) 1 / 2 . The reduced density matrix then 
has off-diagonal elements of the form exp(—c/cr 2 ), which 
vanish as U increases to the Mott limit. For U < 
and \U\N > 2, the function VE 1 ^ can be approximated 
by two gaussians centered at separate points around 
Ni = N/2. As the two peaks sharpen, we converge to the 
two- eigenvalue case described for U — > — oo. Fig.[2]shows 
that S n changes rather sharply around U ~ —2/N, for 
large TV. 

To summarize our findings from the bosonic model, we 
note that the Mott state for U > and Schrodinger cat 
state for U < both possess particle entanglement. We 
have thus identified two generic mechanisms for generat- 
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FIG. 2: (Color online.) 1- and 2- particle entanglement en- 
tropies for TV bosons in two sites. 



ing particle entanglement in itinerant systems. 

Spinless fermions in one dimension — We will now 
consider the ID t-V model; N spinless fermions on L 
sites with periodic boundary conditions: 
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We will use t — 1 units. For repulsive interactions at half 
filling (N = \L), this model has a quantum phase tran- 
sition at V = 2, from a Luttinger- liquid phase at small V 
to a charge density wave (CDW) phase at large V. We 
will focus mainly on V > 0. This model is solvable by 
the Bethe ansatz; however, calculating particle entangle- 
ment entropies S n using the Bethe ansatz is a nontrivial 
problem which we do not address here. 

Limits — For V = (free fermions), the ground state 
is simple in terms of momentum-space modes: a Slater 
determinant of the iV fermions occupying the N lowest- 
energy modes. The n-particle reduced density matrix 
has C(N, n) equal eigenvalues, so that S n = In [C(N, n)], 
independent of the lattice size L. 

In the infinite- V limit, the ground state and hence par- 
ticle entanglement can be simply understood for the case 
of half filling, N = \L. The ground state is an equal 
superposition of two 'crystal' states, and each of them 
gives a separate contribution to the reduced density ma- 
trix. The reduced density matrix has rank 2C(N, n) and 
equal eigenvalues: S n — In [2C(N, n)]. In the notation 
of Eq. ([3]), the subleading term a n intrapolates between 
a n = at V — and a„ — > In 2 for V — > oo for half 
filling. The interpolation details depend on n and N. 

Numerical results — For half-filling (N = ii), Fig. [3] 
presents S n (V), calculated from wavefunctions obtained 
by direct numerical diagonalization. The S n (V) func- 
tion evolves from S FF = In [C(N, n)\ to In [2C(N, n)] ~ 
5ff + 0.6931. For n > 1, we see non-monotonic behavior 
in some cases. At present we have no detailed under- 
standing of the states or particle entanglements at finite 
nonzero V. 

As in our bosonic model, we see Schrodinger cat 
physics in the t-V model also: the V = +00 ground 
state is a superposition of two CDW states of the form 
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FIG. 3: (Color online.) n = 1, n = 2 and n = 3 entan- 
glement entropy in half- filled t-V model (JV = L/2). The 
frec-fermion contribution In [C(N, n)] has been subtracted off. 
The n = 1 plot also displays the effect of a symmetry-breaking 
ec[ci term, with e = 0.1. Inset: position of the maximum as 
function of e. 



|101010...10) and |010101...01). The fragility of this cat 
state can be seen by adding a single-site potential, ecjci, 
or a staggered potential, e' J2i c 2i°2i- The ground state 
then collapses to a single crystal wavefunction, and S n 
drops to In [C(N, n)] (Fig. [3] top panel). 

Phase transition — The particle entanglement entropy 
shows no strong signature of the phase transition at 

V = 2. This is also true after extrapolating to the 
N — > oo limit. (The extrapolated curves are very close 
to the largest- N curves displayed in Fig. [3] and so are 
not shown.) The lack of transition signature is not too 
surprising, because in the definition of the particle en- 
tanglement, the notion of distances (space) enters rather 
weakly. Thus S n is not sensitive to characteristics of 
phase transitions, such as diverging correlation length or 
large-scale fluctuations. 

Away from half-filling — For iV^L/2, the behavior is 
qualitatively similar to the half-filled case, a n increasing 
from zero to an 0(1) value as V increases from zero to 
infinity. (Fig.0J) However, there is no simple picture for 
the V — ► oo limit. Also, a n (V) appears to be monotonic, 
perhaps because a n (V — > oo) is not constrained as in the 
half-filled (CDW) case. 

Note that, except for S n —i in the half-filled case, the 
particle entanglement never saturates the upper bound, 
In [C(L,n)], dictated by Hilbert space size. 

Negative V — An attractive interaction causes the 
fermions to cluster. In the V — > — oo limit, the ground 
state is a superposition (cat state) of L terms, each a clus- 
ter of the N fermions. The cat state can be destroyed as 
in the positive-!^ case. For half-filling with even N, the 

V — > — oo wavefunction yields Si = InN + In 2. There 
are 0(N~ 1 ) corrections for odd N = L/2. 

Eigenvalue spectrum ( majorization ) — The full eigen- 
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FIG. 4: (Color online) (a) S n for N = 7, L = 12 / 
2JV. Horizontal lines are corresponding maximal bounds 
ln[C(L,n)]. (b) negative V, half-filling. Free-fermion con- 
tribution In [C(N,n)] has been subtracted off in each case. 



value spectrum of the reduced density matrices (p n ) of 
course contain more information than the S n alone. For 
n = 1 where S%(V > 0) is monotonic, we numerically 
observe 'majorization' (e.g., Ref. [Io[) of spectra. Obvi- 
ously, there are many other aspects of the full spectra 
that remain unexplored. 

Bosons — An important issue concerns bosonic sys- 
tems which have partial condensation into a single mode, 
so that the leading asymptotic term is not In N. We have 
treated one example, closely related to the fermionic t- 
V model: hard-core bosons on a ID lattice (forbidden 
multiple occupancy, U — oo) with nearest-neighbor in- 
teraction V. The point V — —2 has a 'simple' ground 
state [ll|, which we exploited to find 

S n = vn hxN + O(N ) 

where v = N / L is the filling fraction. A natural interpre- 
tation is that the pre-factor represents the un-condensed 
fraction. Whether this is generic for bosonic systems with 
partial condensation remains an intriguing open ques- 
tion. 

Correlations in subleading term — The canonical re- 
lation S n (N) = In [C(N,n)] + a n allows us to formulate 
correlation effects in terms of the function a n . For free 
fcrmions, for CDW states of the t — V model, and for 
Laughlin states 0, 0| , we have 

a n (FF) = 0, a„(CDW) =ln2, a„(FQH) = n\nm . 

We note that states which are intuitively 'more nontriv- 
ially correlated' have stronger n-dependence in a n . This 
strongly suggests that the a n function is a measure of 
correlations in itinerant fermionic states. It is natural to 
conjecture that the linear behavior of a n is symptomatic 
of intricately correlated states like quantum Hall states, 
and that in generic itinerant states a n will have sublinear 
dependencies on n. 

Equal partitions — In addition to the n <C N behav- 
ior we have focused on here, another promising quantity 
is S n= pj/2- I n Ref. w e presented close bounds for this 
quantity, showing that for fractional quantum Hall states 



of given filling S n —^/2 tends to be higher for more corre- 
lated states. For example S n —jf/2 f° r a Moore- Read state 
is higher than that for a Laughlin state. 

Conclusions — Particle entanglement is an emerging 
important measure of correlations in itinerant many- 
particle quantum systems. In this work, we have set 
the framework for future studies of the asymptotic be- 
havior of particle entanglement. We have also explored 
these quantities in relatively simple itinerant models. We 
have pointed out several different mechanisms for parti- 
cle entanglement in itinerant quantum states, such as 
localization, Schrodinger cat states, and of course anti- 
symmetrization of fermionic systems. Since particle en- 
tanglement is a relatively new quantity on which little in- 
tuition is available, these results will form a much-needed 
basis for future studies. 

Our work opens up a number of questions. Our consid- 
erations have led to an intriguing speculation for bosonic 
systems, relating the leading term in the asymptotic 
(N — > oo) expression for S n (N) to the extent of Bose con- 
densation. A thorough study, addressing several bosonic 
systems, is clearly necessary. In the same asymptotic 
form, one would also like to have a detailed characteriza- 
tion of how the subleading term a n describes correlations. 
More concretely, one could ask "how correlated" a state 
needs to be, in order to have a linear a n function. 
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